Competing Orders and Disorder-induced Insulator to Metal Transition in Manganites 
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Effects of disorder on tlie two competing pliases, i.e., tlie ferromagnetic metal and tfie com- 
mensurate cfiarge/lattice ordered insulator, are studied by Monte Carlo simulation. The disorder 
suppresses the charge/lattice ordering more strongly than the ferromagnetic order, driving the com- 
mensurate insulator to the ferromagnetic metal near the phase boundary in the pure case. Above 
the ferromagnetic transition temperature, on the contrary, the disorder makes the system more in- 
sulating, which might cause an enhanced colossal magnetoresistance as observed in the half-doped 
or Cr-substituted manganites. No indication of the percolation or the cluster formation is found, 
and there remain the charge/lattice fluctuations instead which are enhanced toward the transition 
temperature. 
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The mechanism of the colossal magnetoresistance 
(CMR) in manganites is one of the central issues in 
thephysics of the strongly correlated electronic systems 
B @] ■ By applying the external magnetic field of a few 
Tesla, the resistivity changes of the order of IC^-IO^. Al- 
though the double-exchange (DE) mechanism gives a rea- 
sonable description for the negative magnetoresistance 
[3, it is not enough to explain such a huge response. 
The first proposal for the mechanism of the CMR was 
based on the transition from small to large polaron driven 
by the magnetic field The idea is that the effective 
bandwidth controlled by the DE interaction is enhanced 
by the spin ferromagnetic alignment under the external 
magnetic field, which leads to the crossover from small 
to large Jahn- Teller polaron. However, the change of the 
effective bandwidth is rather small, and this mechanism 
does not lead to such a huge effect. Thus, this single- 
particle picture is not satisfactory to explain the CMR. 

Instead, recently, it has been recognized that the CMR 
is a collective phenomenon due to the many-body cor- 
relation. One of the present authors proposed a new 
mechanism of CMR based on the multicritical fluctuation 
between the ferromagnetic metal (FM) and the charge- 
ordered insulator (COI) occurring near x = 0.5 (x: hole 
concentration) There, the enhanced fluctuations near 
the multicritical point trigger a giant response to the ex- 
ternal magnetic field. An evidence for this multicritical 
scenario is the scaling law for the magnetization curve 
(Arrott plot). However, this scaling analysis holds only 
for some class of materials, while the others do not show 
a good agreement. 

In the latter class of materials, the disorder appears 
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FIG. 1: Phase diagram of ini/aBai/aMnOa from Ref. [8]. Tc, 
Icq, 7n, and Tsg denote the ferromagnetic, charge-ordering, 
antiferromagnetic, and spin-glass transition temperatures, re- 
spectively. The subscript '(order)' are for the (in,Ba)-ordered 
materials. The shaded area shows the charge-ordered insulat- 
ing state in the ordered materials. The hatched area repre- 
sents the disorder-induced ferromagnetic phase. 



to play an important role. Especially the spatially inho- 
mogeneous structure has been observed experimentally 
and the percolation mechanism of the CMR is phe- 
nomenologically proposed based on this observation 0. 
Namely it is assumed that the resistivity is determined by 
the percolating conduction paths, which are so sensitive 
to the pattern of the coexisting metallic and insulating 
regions under the influence of the disorder. The external 
magnetic field changes this pattern, and hence causes a 
large change of the resistivity. 

Recent experiments in the half-doped manganites 
LTii/2Bai/2Mn03, where Ln is a rare-earth element, have 
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elucidated more explicitly the importance of the disorder 
d, . In these compounds, the strength of disorder due 
to the solid solution of (Ln,Ba) ions can be tuned by 
careful treatment of the synthesis |0| . Figure ^ shows 
the experimental phase diagram reported in Ref. Q. In 
the 'clean' limit where Ln and Ba ions form a periodic 
layered structure, the phase diagram shows a multicrit- 
ical behavior where FM and COI compete with each 
other. On the other hand, in the 'disordered' case where 
the arrangement of (Ln,Ba) ions is random, the phase 
diagram changes in a very asymmetrical manner: FM 
is suppressed but survives at finite temperatures, while 
COI disappears and instead some glassy state is realized 
at low temperatures. This leads to a nontrivial regime 
where the disorder induces the transition from COI to 
FM (hatched area in Fig.P). The enhanced CMR effect 
is observed above Tc in this regime |8j| . Similar phenom- 
ena are observed also in Lno.55(Ca,Sr)o.45Mn03 and 
in Lno.sCao.sMnOs with the Cr substitution of Mn T^. 

Therefore, to understand the mechanism of the CMR, 
it is crucial to clarify the effects of the quenched disor- 
der on the competition between the FM and CO, and 
on the metal-insulator transition, which we undertake in 
this letter. Employing the unbiased numerical method 
to treat large fluctuations properly, we reveal that these 
two orders show contrastive responses to the random po- 
tential which indeed leads to the disorder-induced tran- 
sition from COI to FM. Through the systematic study of 
density of states, conductivity, and fluctuations, we elu- 
cidate the mechanism of the enhanced CMR effect in this 
regime. 

We consider a minimal model to describe a competition 
between FM of the DE origin and COI concomitant with 
the lattice distortion, whose Hamiltonian reads 

{i])<j i 
I i (tj) i 

Here, the first two terms represent a simplified DE model 
which consists of the electrons' nearest-neighbor hopping 
t and the Hund's-rule coupling Jh . We consider the Ising- 
type Hund's-rule coupling in the limit of Jh — ^ oo for 
simplicity, which retains an essential physics of the DE 
mechanism The third term describes the electron- 
phonon coupling. Here, we consider the classical phonons 
which couple to the electron density (breathing mode). 
The next two terms are the elastic energy of phonons. 
We take account of the cooperative aspect of the lattice 
distortion by the fifth term which reflects the fact that a 
displacement of an oxygen affects both MnOg tetrahedra 
connected by the oxygen. The last term incorporates the 
random on-site potential energy £i which couples to the 
electron density rii = The alloying effect of A- 

site ions in j4Mn03 and the Cr substitution into the Mn 



site cause the structural and electrostatic disorder, which 
disturbs the on-site potential through the Madelung en- 
ergy. Since large fluctuations in the competition between 
different orders are crucial in this study, we employ the 
Monte Carlo (MC) method in which the configurations 
of {Si} and {Qi} are stochastically sampled In the 
presence of the disorder, we take the random average for 
different configurations of {si}. 

We consider model ^ on the 2D square lattice in the 
half-doped case with x ~ 0.5 (0.5 electron per site on 
average). We set the half-bandwidth W = At as an en- 
ergy unit, and take A = 0.1 throughout this paper. A 
similar model was studied in a different context for the 
low- a; regime in the absence of the disorder . We con- 
sider the binary-type distribution of the random poten- 
tial, Ei — ±A. We apply the systematic analysis on the 
finite-size effect by using the series of L x L-site systems 
with L = 4, 6, 8 to distinguish the long-range ordering 
from short-range correlations. The numerical details will 
be reported elsewhere [l^ . 

Figure El (a) shows the phase diagram of model 
in the absence of the disorder (A — 0). Note that the 
abscissa is reversed for comparison with experimental re- 
sults. At g = 0, the electronic part of model ^ becomes 
a pure DE model which exhibits the FM state below 
T(^"' (the superscript (0) represents the case of A = 0). 
As g increases, T^^ decreases and the checkerboard-type 
COI concomitant with the cooperative lattice distortion 
appears below T^q. In the intermediate region where 



Tq"'' and Tqq intersect, we have the phase with both the 
ferromagnetism and the charge/lattice ordering (F-f-CO 
phase). Thus, the phase diagram exhibits the tetracrit- 
ical behavior where four different phases meet at one 
point. The density of states (DOS) shown in Fig. [3 (d) 
indicates that the phases below T^q including the F+CO 
phase have a finite energy gap. This tetracritical behav- 
ior can be turned into the bicritical one when we include 
a competing term between the ferromagnetism and the 
charge/lattice ordering by hand 'iG^I, and the following 
results on the disorder effects are qualitatively similar in 
these multicritical phenomena. 

The emergence of the F-I-CO phase is due to the ro- 
bustness of the DE ferromagnetic interaction even in the 
insulating state. The DE interaction is effective if the 
electron wavefunction spreads, at least, for a finite spa- 
tial range. We note that this phase is al so p redicted in 
more realistic DE model for manganit es IITI H^. and is 
indeed observed in recent experiments jl9l | . 

When we switch on the random potential, the phase 
diagram shows a drastic change. By introducing weak 
disorder with A = 0.1 (Fig. |21 (b)), Tco is suppressed 
remarkably, whereas Tq is lowered only slightly. Accord- 
ingly, the tetracritical point shifts toward larger g and 
lower T. For stronger disorder with A = 0.2 (Fig. [21 
(c)), Tqo disappears down to the lowest temperatures in 
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FIG. 2: (a)-(c) Phase diagram of model (^3 for different 
strength of the random potential (a) A = 0.0, (b) A = 0.1, 
and (c) A = 0.2. The lines are guides for the eyes. The 
gray data and lines in (b) and (c) are the results in (a) for 
comparison, (d)-(g) Density of states at T = 0.016 for L = 8 
clusters: (d) A = 0.0 for various g; (e) g = 0.4, (f) g = 1.6, 
and (g) g = 0.8 by varying A. The Fermi energy is at lj = 0. 
Only the typical error is shown for each case. 



the present MC calculations (T ;^ 0.008). This clearly 
shows the fragility of COI against the random potential 
compared to the DE ferromagnetism 0| . As a result of 
this asymmetric effect, there appears the regime where 
COI is turned into FM by introducing the random po- 
tential, which is not clearly found in the previous study 
[l7|. These features are in good agreement with the ex- 
perimental results 0, 0| . 

The DOS gives insights on the nature of these disor- 
der effects. As shown in Fig. |21(e) and (f), DOS in both 
the small-g FM regime and the large-g COI regime are 
not much affected by the random potential. Especially, 
Fig. 121 (f) indicates that the energy gap in the large-g 
regime persists even when the charge/lattice long-range 
order is destroyed by the disorder. This means that a 
short-range correlation remains and the local lattice dis- 
tortion persists to open the gap. Therefore, the collapse 
of the charge/lattice ordering is not due to the disappear- 
ance of the amplitude of the lattice distortion but due to 
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FIG. 3; Optical conductivity for g = 0.8 at (a) T = 0.02 and 
(b) T = 0.06 (L = 6). (c) Average of the square of the lattice 
distortion {L = 8). See the text for details, (d) Schematic 
picture of the temperature dependence of the resistivity. 



the disturbance of the phase of the commensurate order- 
ing with the wavenumber (tt, tt). The random potential 
acts as the 'random field' to the CO order parameter, 
and its phase coherence is lost. On the contrary, the fer- 
romagnetic order is robust against such pinning since Tc 
is essentially determined by the kinetic energy of elec- 
trons. These give a comprehensive understanding of the 
asymmetrical effect of the disorder on the phase diagram. 

The most remarkable behavior of DOS is obtained near 
the phase boundary between the FM state and the F+CO 
insulating state. In the F-I-CO phase but close to the 
phase boundary at A = 0, as shown in Fig.|2](g), a small 
but clear gap in DOS below T^q is rapidly collapsed by 
introducing the disorder and the system appears to be 
metallic. This indicates that the random potential easily 
destroys COI and induces FM in the proximity of the 
metal-insulator transition. This result also agrees well 
with experimental results which strongly suggests 

that the simple model captures the essential physics 
of this competition between FM and CO in manganites. 

Above Tc of the disorder-induced metaUic state, the 
CMR is found to be much enhanced experimentally 
There, the system is insulating and shows a large 
drop of the resistivity near T = Tq toward the low-T 
metallic phase. Our results are satisfactorily consistent 
with this behavior: FiguresEKa) and (b) show the optical 
conductivity a{uj) for T < Tc and T > Tq, respectively, 
which are obtained by the Kubo formula. In the low-T 
regime, a{uj) shows a disorder-induced insulator to metal 
transition which is consistent with DOS in Fig. |3 (g). 
On the contrary, at the high-T regime, the disorder sub- 
stantially reduces the low-energy weight and develops a 
quasi-gap structure in cr{uj), namely, it makes the system 
more insulating. These tendencies are also indicated in 
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FIG. 4: (a) Temperature dependence of the standard devi- 
ation of the square of the local lattice distortion at A = 0.2. 
Transition temperatures for FM (COI) in the absence of the 
disorder are indicated by the gray (black) arrows, (b) Sus- 
ceptibility for the (7r,7r) lattice ordering at g = 1.0. Both are 
for L = 8 clusters. 

the magnitude of the lattice distortion, (X^i )/-^^' 
shown in Fig.|31(c), which correlates with the localization 
of electrons. The distortion increases for T > Tq while 
it decreases for T < Tq when the random potential is 
introduced. These contrastive effects of the disorder sug- 
gest the temperature dependence of the resistivity p as 
schematically drawn in Fig. 13 (d) . In the absence of the 
disorder, the system causes the COI transition at T(1q, 
where p shows a sharp upturn. The disorder enhances the 
insulating nature in the high-T regime, hence p becomes 
larger while it does not show any anomaly corresponding 
to the COI transition. Below Tc, the disorder induces 
FM, therefore p should show a large drop at T ^ Tc as 
shown in the figure. The large change of p potentially 
leads to a huge response to the external magnetic field as 
is indeed observed in the CMR manganites (sj, [ij . Thus, 
the origin of the enhanced CMR can be attributed to the 
contrastive influence of the disorder in low-T and high-T 
regimes near the metal-insulator phase boundary. 

We discuss here this insulating state for T > Tc in 
more detail. In the percolation scenario for the CMR 
0, the system is supposed to be a coexisting mixture of 
the metallic and the insulating clusters below the FM or 
COI transition temperatures in the absence of the disor- 
der, and T(^q, respectively. However, our simulation 
does not show any clear indication of this coexistence. 
Figure ^ (a) shows the standard deviation of the square 
of the lattice distortion {Ql) among the sites as a func- 
tion of the temperature. The standard deviation changes 
monotonously and remains small (compare to the values 
in Fig. 131(c)). This is contrary to the percolation picture 
in which we expect an enhancement for T < T^'^ or T^q 
due to the formation of the metallic (smaller (Qf )) and 
the insulating (larger {Qf)) clusters. Instead of the clus- 
ter formation, we find fluctuations or short-range corre- 
lations of the charge/lattice ordering. Figure^ (b) shows 
the susceptibility for the (tt, tt) lattice ordering. In the 
absence of the disorder, the susceptibility shows a sharp 
cusp at T = T^Q, which leads to the divergence in the 



thermodynamic limit T ^ oo. Although the disorder 
destroys COI and makes the cusp obscure, the fluctua- 
tion remains finite and is enhanced as T Tq. The 
contrastive behavior between above and below Tc in the 
disordered case is in accord with this temperature de- 
pendence of X- Namely the CO fluctuation is enhanced 
towards Tc from above while it is suppressed at the low- 
est temperature. This could be regarded as the reminis- 
cence of the multicritical phenomenon in the pure case. 
This is also in good agreement with the diffuse scatter- 
ing experiments which indicate a l arg e fluctuation of the 
charge/lattice ordering above Tc [l^illyl- We note that 
this fluctuation is suppressed as A increases, which may 
account for the deviation from the scaling law ^] . There- 
fore, our results indicate that the phase above Tc is not 
likely the inhomogeneous mixture of the static clusters 
but a rather homogenous state with the thermodynamic 
charge/lattice fluctuations. 
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